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Abstract 

We investigate the correlation between nuclear neutron radii and the radius of neutron stars. 
We use a well-established hadronic SU(3) model based on chiral symmetry that naturally includes 
non- linear vector meson and scalar meson - vector meson couplings. The relative strengths of the 
couplings modify the nuclear isospin-dependent interactions. We study the dependence of nuclear 
and neutron star radii on the coupling strengths. The relevance of the results for parity- violating 
electron- nucleus scattering and the URCA process in neutron stars is discussed. 
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I. INTRODUCTION 



In light of the new and planned radioactive beam facilities the study of neutron-rich 
nuclei has received increased attention in recent years. Completing the chart of metastable 
nuclei towards the proton and neutron drip lines is the major goal of the endeavor with 
important information for the stellar nucleosynthesis. From a theoretical point of view 
nuclear calculations toward the neutron drip-line are still quite problematic in part due to 
the rather poor knowledge of the isospin-dependence of the nuclear forces. Without going 
to extremely neutron-rich nuclei there is another complementary information by studying 
nuclear properties and the ultimate neutron-rich systems, i.e. neutron stars, in a combined 
approach. The radius of a neutron star is sensitive to isospin forces and over a significant 
part of the star the corresponding nuclear densities have values that also occur in atomic 
nuclei. 

A rather sensitive measure for isospin effects in a normal nucleus is its neutron radius 
compared to the proton radius. In an upcoming experiment at Jefferson Laboratory an ac- 
curate measurement of the neutron radius of 208 Pb is planned using parity-violating electron 
scattering 

As a basis for studying correlations between the nuclear radii with neutron star radii a 
theoretical approach has to be employed that covers nuclear structure as well as neutron 
star calculations in a unified way. Since the properties of neutron stars can change signif- 
icantly by including hyperon degrees of freedom, a SU(3)-flavor approach is a convenient 
choice for such an investigation. In order to study the dependence of nuclear and neutron 
star radii we investigate the variation of the radii by modifying terms that influence the 
isospin dependence and nuclear skin effects. Here we follow the analysis of |2j, ||[ where in 
a Walecka-type model the sensitivity of nuclear and star radii on the coupling strengths 
between isovector vector mesons and scalar and vector mesons was investigated. As it turns 
out such an analysis can be done in our approach in a very natural way. 

In the following we present a calculation of the radii of neutron stars and the 208 Pb proton 
and neutron radii in a framework that is based on a model constructed from chiral SUp(3) 
arguments. The outline of the article is as follows. First in Section II we briefly introduce 
the general framework of the chiral model looking at several meson interaction terms. Then 
in section III we discuss the results of the calculation, studying the dependence of neutron 
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star radii and nuclear neutron skins on the various coupling constants. From those results we 
derive the direct correlation between neutron skin and neutron star radius. As a test of our 
model we also look at the proton skin of proton-rich argon isotopes. Finally, we investigate 
the consequence for neutron star cooling via the URCA process. 

II. THEORETICAL FRAMEWORK 

In our nuclear structure as well as in the neutron star calculation we adopt a generalized 
flavor-SU(3) a — u type model. The reason for adopting this model is two-fold. For one 
it includes, by construction, a nonlinear coupling of the p and uj meson as well as coupling 
between the vector mesons and scalar mesons, so there is no need to add additional terms 
to an established model to study the type of isospin dependencies mentioned before. Sec- 
ondly, in addition to successful applications in one- and two-dimensional nuclear structure 
calculations [|], ||, the same model also shows a sensible behavior at high temperatures and 
densities||, which is quite important in this discussion as the inclusion of nonlinear terms 
in a calculations involving matter in extreme regions of temperature or density, has the gen- 
eral problem that this can produce rather uncontrolled extrapolations. A phenomenological 
check of the soundness of the model predictions at extreme values is therefore reassuring. 
The model used here has been successfully tested in all those regimes. 

The model includes the complete set of the lowest SU(3) multiplets of baryons and mesons. 
In addition extensions to higher multiplets have been discussed 0. 

Restricting ourselves to those degrees of freedom relevant for the discussion in this paper 
the general SU(3) Lagrangian || reduces to the following structure: The interaction term 
of baryons, mesons, and the photon reads 



where now the SU(3) baryon octet B is reduced to the isospinor (£). The various coupling 
constants of mesons and baryons result from the SU(3) structure || and the interaction of 
the scalar fields follows as 
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Here the various scalar fields a, 5, and ( correspond to the non-strange isoscalar (a), isovector 
(5) quark-antiquark states, and the strange anti-strange state (£). The field % is the scalar, 
isoscalar glueball field x introduced in || [HJ. xoi a Q an d Co (^o — 0) denote the vacuum 
expectation values of the fields, which are generated via spontaneous symmetry breaking. 
The term 

£esb = - (j^j [xa + y(] (3) 

introduces explicit chiral SU(3) symmetry breaking. The coupling strengths are x = Tn 2 f n 
and y = y/2m^f K — ■^=rri^f 7r . f n and fx are the pion and kaon decay constant, respectively. 
The vector meson self-interaction terms have the structure 

Cvec,NL = -a (Tr [V^]f - 6Tr [(V^) 2 ] (4) 

where V denotes the 3x3 matrix of the vector meson multiplet. In the case of nuclei this 
reduces to the form 

C vec = -g\ (cu A + + P 4 ) + • • • ( 5 ) 

with the zeroth components of the isoscalar u and neutral isovector p field. The couplings 
g\ and (3 are related to a, b of equation [| via 

In addition the vector mesons interact with the glueball field x an d the nonet of scalar 
mesons S via 

C vs = c X 2 Tr[^y M ] + dTr[£V%] (7) 

which reduces to 

C vs = [c X 2 + 2d(a 2 + S 2 )] (lu 2 + p 2 ) + . . . (8) 

In the numerical calculation we investigate the relative strength of the two couplings of the 
vector mesons to the glueball field as well as to the scalars by varying the quantity r a defined 

as 



Of. 



(1 _ r ^ 5 + rXl (l-^ + r^ + ^(^ + /) (9) 

where the factor a was determined in the general parameter fit xm with r a = H. r a 
shifts the relative importance of the gluonic and quark-type scalar fields for generating the 
masses of the vector mesons. The mass of the baryons are dynamically generated via their 
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interaction with the scalar fields following Eq. (|]), which automatically leads to density and 
temperature-dependent masses. Explicitly, the baryon masses m* are given by 



9ia<? + 9is5 + gi(C ■ (10) 



In general there is no unique way of introducing nonlinear vector couplings in a hadronic 
model. Here we introduce the simplest structures within our framework generating the 
interaction terms, in fact those terms had been introduced well before this specific investi- 
gation of isospin effects started. The non-linear couplings in the u channel turn out to be 
very important to obtain a reasonably soft equation of state of symmetric nuclear matter, 
leading to acceptable value for the compressibility of k < 250 MeV. In addition, away from 
the saturation density the equation of state of neutron matter is soft enough to generate 
neutron stars with masses M < 2M so i ar [11]. 



III. RESULTS 

A. Neutron skins 

Here we show the numerical results obtained by solving the equations of motions derived 
from the Lagrangian terms discussed above in a mean field approximation. In particular we 
will vary the two coupling strengths /3, Eq. (||), and r a , Eq. (||) and study their influence 
on various observables. 

First we determine the rms radii of neutrons and protons in 208 Pb. In addition as com- 
parison we also look at the case of the singly magic 138 Ba which might be of interest in 
future parity-violation experiments. Barium isotopes have been suggested as candidates for 
atomic parity violation experiments |T3 in analogy to the cesium experiments ||14j|. In those 



experiments uncertain neutron distributions are a major part of the theoretical error. The 
calculation of 138 Ba is supposed to give an idea of the range of uncertainty for neutron skins 
in barium. 

The equation of motions are solved assuming spherical symmetry using the established 
parameter set \m fU- The calculation is repeated by changing the two coupling parameter 
P and r a . The resulting values for the difference of proton and neutron radii Ar np = r n — r p 
are shown in Figure In the variation of the parameters we follow two general procedures. 
First we use the parameters \m an d only modify the coupling strength (3 (which has a 




FIG. 1: Difference of neutron and proton radius in 208 Pb as function of nonlinear vector coupling 
strength (3 with (squares) and without (circles) refitting of parameters (r a = 0). 

value (3 = 1 in Xm), without simultaneously modifying any other parameters. In a second 
run we adjust the other parameters that mainly influence the isospin behavior, gf and 
the nucleon-rho coupling g^ p to ensure a correct description of the lead charge radius and 
binding energy. The result of both approaches is shown in the figure. The general behavior 
is the same for both approaches. The best-fit estimate for neutron-proton radius difference is 
Ar np = 0.25 fm. This number is in general agreement with other models based on relativistic 
descriptions of nuclear matter [IS] and is bigger than typical Skyrme-based non-relativistic 



projections of Ar ra „ < 0.2 fm [fig] . There is a significant dependence of Ar np on the coupling 



strength while varying (3 between about — 1 < (3 < +1.4, which is about the maximum range 
in which one can obtain stable nuclear solutions. We will look at the parameter range in more 
detail. In order to obtain a more believable quantitative result a refitting of the parameters 
is necessary, reducing Ar np by up to 20 percent. In the following, if not noted otherwise, 
we will always perform a readjustment of the parameters as explained above. Figure |2| 
shows results for Ar np varying the coupling r a for different values of (3. The shaded area 
reflects the expected accuracy of the neutron radius measurement of about ±1%. Note the 
wiggle for the line with (3 = 1.4, signalling the onset of proton and neutron rearrangements 



6 




FIG. 2: Difference of neutron and proton radius Ar np in Pb as function of relative vector- 
scalar meson coupling strength r a for different values of (3 = {—0.6,-0.2,0.2,0.6,1.0,1.2,1.4}. 
The shaded area denotes the expected experimental accuracy for the r n measurement (tentatively) 
centered at Ar np = 0.25 fm. 

in the solutions of the nuclear equations that violate experimentally known nuclear charge 
distributions. 

We repeat the same calculation for the case of 138 Ba, Fig. |3]. The results are essentially 
similar. The dependence of the neutron skin is somewhat less pronounced compared to the 
lead case. 

B. Nuclear matter 

In order to check that aside from a reasonable description of 208 Pb the range of parameter 
values yields adequate behavior of nuclear matter at saturation we study the values for the 
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FIG. 3: Difference of neutron and proton radius in 138 Ba as function of coupling strength r a . 
Curves for different values of (3 = {—0.6, —0.2,0.2,0.6, 1.0, 1.2} are shown (cf. Fig. §). 

compressibility k and the asymmetry energy 0,4 at saturation density po, which has a value 
p ~ -15 GeV/fm 3 for all studied parameter sets, k and a 4 are defined as 



- ^0 p2 \p=po, 04 " 2 ^ _ pp)2 Ip=po ■ ) 



Fig. ^ shows the r CT dependence of k. The two curves correspond to two extreme values 
of P, the other results lie between these curves. Given these numbers one cannot exclude 
any parameters within the range under consideration, all give reasonable values for the 
compressibility of nuclear matter between about 205 and 225 MeV. The results for the 
asymmetry energy are more restrictive. Assuming that acceptable value for the asymmetry 
should be somewhere between 28 and 38 MeV, values of r a larger than 0.2 yield asymmetries 
04 that are too small, whereas given small couplings of scalar and vector mesons the coupling 
term (3 should be larger than 0. However, a quite large range of reasonable parameter values 
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FIG. 4: Results for the compressibility k as function of the coupling strength r a . The two curves 
show the results for the different values of (3 = 1.2 (full line) and j3 = —0.6 (dotted line). 

persists. 

C. Neutron stars 



As a complementary source of information on the isospin-dependence of nuclear forces a 
neutron star probes the extreme neutron-rich system with a content of about 90% neutrons 
of all baryons in the star. We have performed calculations of neutron star properties in this 



model before in 11, 12]. The calculation is done by integrating the Tolman-Oppenheimer- 



Volkov (TOV) equations [lTj for the star given an equation of state e(P) generated from the 
model Lagrangian: 

dP(r) 



dr 



(e+P)(47rr 3 P+m) 
r 2 (l-2m(r)/r) 



(12) 
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FIG. 5: Results for the asymmetry energy 04 of saturated nuclear matter as function of the 
coupling strength r a . The horizontal axis shows the value for r a whereas the different curves are 
distinguished by their values of (3 ranging from —0.6 to 1.2 (see Fig. |||). 

^P- =47rrMr) , (13) 
dr 

where P and e are the pressure and energy density of the nuclear matter. m(r) is the 
gravitational mass of the star inside of the radius r. By fixing the initial value for the 
pressure P c in the center of the star, with the knowledge of the equation of state e c (P c ) one 
can integrate equations (0) and (|13|) from the center of the star outward to the point of 
vanishing pressure, which defines the star radius R. The calculation of e(P) is discussed in 
detail in |J. The outer crust of the neutron star is modelled as discussed in []TTj] . 

Figure ^ shows the resulting neutron star radii for a typical neutron star with a mass 
of 1.4 solar masses. As in the nuclear case we checked for the consequences of refitting 
the parameters. But in this case the refitting does not affect the neutron star properties 
appreciably as can be inferred from the figure. There is a strong dependence of the neutron 
star radius on the coupling strength with a variation of about 5.5 km over the whole range 
of values for f3 that yield reasonable nuclei. However, with reference to the results for the 
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asymmetry shown in Fig. [5] in this case for r a = (3 should be chosen to be larger 0, which 
limits the variation of the star radius to about 1 km. Another uncertainty entering the 
theoretical determination of neutron star radii comes from the treatment of the hyperons in 
the star. Although our model is based on a SU(3) model and therefore automatically contains 
hyperons (A, £ and H baryons, see [O] for a detailed discussion), due to the poor knowledge 



of the behavior of hyperons in dense matter, results might vary with the treatment of the 
hyperon sector. To check the limit of the uncertainties we switch off (by hand) the hyperon 
densities and repeat the radius calculations for the example of the best-fit parameter set 
Xm- The variation for a 1.4 solar mass star, Rh ype r = 11-42 km and R no hyper = 12.46 km 
is of the same order as the variation due to the meson-meson couplings. So in addition to 
a good knowledge of the nuclear isospin forces some information on hyperons in matter is 
needed to draw more reliable conclusions from neutron star radius data. 

A complete set of correlations between Ar np and R ns is given in Fig. [F[ The various 
lines are results for different choices of (3 and show the dependence on the scalar-vector 
meson coupling r a . The result shows that the information contained in (3 and r a is rather 
complementary. Whereas the star radius largely depends on (3 and little on r a in the case 
of the neutron skin the result is the opposite, a strong r a dependence and a small influence 
from the (3 value. However, as can be seen from the hatched area a 1% measurement of the 
neutron radius of lead covers about the whole range of results except for very small (3 and 
large r a values which are largely excluded by the results for the asymmetry energy. 

D. Additional tests 

As a cross check to our discussion we looked at the proton skin in different isotopes 
of argon. The skin was experimentally deduced by combining measurements of heavy-ion 
reaction cross sections|18[ with atomic isotope shift data [ I9| . The resulting values, with 



rather large error bars, are shown in Figure |8], compared to the theoretical calculations using 
the same model for several values of the nonlinear coupling varying (3 (r a is set to 0). One 
can see that with the exception of the point at the high-end of coupling strengths for the 
32 Ar isotope the agreement with data is quite good, independent of the f3 parameter. This 
gives some confidence in the analysis that a change of the coupling does not automatically 
lead to distorted results of other isospin-related observables. 
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FIG. 6: Radius of a neutron star with M = 1.4 M so i ar varying the nonlinear coupling strength (5 
for a fixed r a = 0. The solid (dashed) line shows the results with (without) refitted parameters 
(see text). 

Finally, in []2D| the importance of the URCA process for neutron star cooling was discussed 
in this context. The URCA process is thought to be an important mechanism for cooling 
the neutron star after the supernova explosion at reasonably fast rates in agreement with 
astronomical observations. This is done via the process n — > p + e~ + v e where the neutrino 
escapes and carries energy away from the system. This process can proceed when the fermi 
momenta of the particles involved match : 

Pn<Pp+Pe (14) 

Thus for the URCA process to take place one needs a significant amount of protons of 
roughly 11% in the system. As the proton to neutron ratio in the star is influenced by 
the isospin-dependence of the nuclear forces one has to check whether changing the force 
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FIG. 7: Relation of the difference of neutron and proton radius Ar np of 2U8 Pb (full symbols) and 
the neutron star radius (M = 1.4 M so i ar ). The coupling strengths j3 and r a are varied. Results for 
different values of (5 are marked with different symbols with values of 1.2, 1.0 0.6, 0.2, -0.2, -0.6 
from the left to the right. For every value of (5 we varied r a between and 0.3, with large values 
of r a yielding small values or Ar np as can also be seen in Fig. ||. 

might switch off the URCA process. We again consider a 1.4 solar mass star and compare 
the central density of the star with the minimum density for the URCA process to occur 
according to equation 14. We see from Figure ^]that the minimum density is always reached. 



Thus the discussion of the correlation of the star radius and neutron skin is not affected by 
additional constraints from neutron star cooling considerations. 

IV. CONCLUSION 

We studied the influence of couplings between vector and scalar mesons as well as non- 
linear vector-meson interactions and their influence on isospin-related observables following 



the discussion in [|[ (see also [21]). As basic model for our calculation we used a chiral 



a — oj type SU(3) model. This model has the benefits of including those nonlinearities from 
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proton skin in Ar isotopes 
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FIG. 8: Proton skin for various Argon isotopes in comparison to experimental values |]T^]. Results 
for different nonlinear coupling strengths are shown, the squares, circles, diamonds and stars mark 
results for (5 = 1.4, 1.0, —0.2, —0.6, respectively. 

the start and implements them in a rather general SU(3) coupling scheme. Our model has 
been successfully employed in nuclear structure calculations and in calculations involving 
high-density and high-temperature environments. Since we also study neutron stars the 
benefit of having an SU(3) model is that hyperons are automatically included in the model 
which is quite important for getting reasonably small neutron star masses in accordance 
with observation. The results show clear dependence of the neutron skin in lead and barium 
on couplings of vector and scalar mesons, whereas the neutron star radii depend largely on 
non-linear terms in the vector-meson channel. There is a distinct correlation of neutron 
star radii and neutron distributions in nuclei as has already been observed in other nuclear 
structure models @, [|. However, although the upcoming measurements of the lead radius 
using parity- violating electron scattering experiments certainly will give additional informa- 
tion for modelling the isospin-sector of the nuclear forces, it is likely not stringent enough 
to nail down the parameters. Also possible measurements of neutron star radii can only be 
useful in this regard if the accuracy of the radius measurement is well below 1 km. At this 
level of radius uncertainty also the detailed treatment of the hyperons in the neutron star 
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FIG. 9: Critical density for the onset of the URCA process for different nonlinear coupling strengths 
(3 setting r a = 0. The squares show the maximum density for a neutron star with a typical mass 
of M = 1.4 M ' S oiar- The circles represent the minimum density necessary for the URCA process. 

enter, which will make it very difficult to disentangle the different effects on the star radius. 

To complete our investigation we have shown that the requirement of fast neutron star 
cooling does not constrain this discussion. Also the measured value of the proton skins of 
argon isotopes can be reproduced and are not strongly affected by varying the nonlinear 
couplings. 



Acknowledgments 

This work was supported by the U.S. Department of Energy, Nuclear Physics Division 
(Contract No. W-31-109-Eng-38). 



[1] JLab experiment E-00-003, R. Michaels, P.A. Souder, CM. Urciuoli (spokespersons) 



15 



[2] C. J. Horowitz and J. Piekarewicz, Phys.Rev. C64, 062802 (2001). 
[3] C. J. Horowitz and J. Piekarewicz, Phys.Rev.Lett. 86, 5647 (2001). 
[4] S. Schramm, preprint nucl-th/0207060| , Phys. Rev. C. to be published. 



[5] Ch. Beckmann, P. Papazoglou, D. Zschiesche, S. Schramm, H. Stocker, and W. Greiner, Phys. 

Rev. C 65, 024301 (2002) 
[6] D. Zschiesche, P. Papazoglou, S. Schramm, J. Schaffher-Bielich, H. Stoecker, and W. Greiner, 

Phys.Rev. C63 (2001) 025211. 



[7] D. Zschiesche, S. Schramm, J. Schaffher-Bielich, H. Stoecker, W. Greiner, preprint qucL 



| th/020902l , to be published in Phys. Lett. B. 
[8] P. Papazoglou, D. Zschiesche, S. Schramm, J. Schaffher-Bielich, H. Stocker, W. Greiner, Phys. 

Rev. C59, 411 (1998). 
[9] J. Schechter, Phys.Rev. D21, 3393 (1980). 
[10] P. Papazoglou, S. Schramm, J. Schaffher-Bielich, H. Stcker, and W. Greiner, Phys.Rev. C57, 
2576 (1998). 

[11] M. Hanauske, D. Zschiesche, S. Pal, S. Schramm, H. Stocker, and W. Greiner, Ap. J. 537, 958 
(2000). 

[12] S. Schramm and D. Zschiesche, preprint nucl-th0204075, submitted to J. Phys. G. 

[13] M. J. Ramsey-Musolf, Phys. Rev. C60 (1999) 015501. 

[14] S. C. Bennett and C. E. Wieman, Phys. Rev. Lett. 82, 2484 (1999). 

[15] D. Vretanar, G. A. Lalazissis, and P. Ring, Phys. Rev. C 62, 045502 (2000). 

[16] K. Pomorski, P. Ring, G.A. Lalazissis, A. Baran, Z. Lojewski, B. Nerlo-Pomorska, and M. 

Warda, Nucl.Phys. A624, 349 (1997). 
[17] R.C. Tolman, Phys. Rev. 55 (1939) 364; 

J.R. Oppenheimer and G.M. Volkoff, Phys. Rev. 55 (1939) 374. 
[18] A. Ozawa et al., preprint RIKEN-AF-NP-424, May 2002. 
[19] A. Klein et al., Nucl. Phys. A607, 1 (1996). 

[20] C. J. Horowitz and J. Piekarewicz, preprint |nucl-th/020706^ (2002). 
[21] A. R. Bodmer and Q. N. Usmani, Argonne preprint (2002). 



16 



